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We report on the development of a new theoretical tool that allows for isospin projection 
of Slater determinants and we present its first applications. In particular, we determine 
the isospin mixing in ground states of N = Z nuclei and discuss its dependence on 
the size of the harmonic-oscillator basis used in the calculations. We also discuss the 
unphysical contribution to the isospin mixing caused by the spurious isospin-symmetry 
breaking inherent to the mean-field approach. We show that these contributions may be 
as large as 30% of the value of the isospin-mixing parameter. 



1. Introduction 

Self-consistent mean-field (MF) approach is practically the only formalism allowing 
for large-scale no-core computations in heavy open-shell nuclei with many valence 
particles. Inherent to the MF approach is the mechanism of spontaneous symmetry 
breaking, which is essentially the only way to allow for incorporating a signifi- 
cant part of many-body correlations into a single intrinsic (symmetry-breaking) 
Slater determinant. However, unlike in the cases of rotational or translational 
symmetry-breaking schemes, violation of the isobaric symmetry has two distinc- 
tively different sources. The unwanted or unphysical source pertains directly to the 
MF approximation ^ «J ^. It manifests itself very clearly in the ground-state wave 
functions calculated by using isospin-invariant interactions, like Skyrme or Gogny 
forces, with Coulomb force neglected. Indeed, such calculations manifestly break the 
isospin symmetry in all but the N = Z systems, simply because the self-consistent 
proton and neutron wave functions are then different. The second source of the 
isospin-symmetry violation is of strictly physical nature and is caused mostly by 
the Coulomb field and, to a much lesser extent, by strong-force isospin-non-invariant 
components. 

Hereby, we report on the development of a new theoretical tool that allows for 
isospin projection (after variation) of Slate r determinants. It has been implemented 
within the Hartree-Fock code HFODE0^. First, the Slater determinants are deter- 
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mined in a standard way by minimizing the Skyrme functional plus the Coulomb 
energy. Both direct and exchange terms of the Coulomb energy are calculated ex- 
actly. We allow for arbitrary spatial deformations of these intrinsic states. Second, 
the isospin-projected components are determined, and third, they are mixed so as 
to rediagonalize the total Skyrme-plus-Coulomb Hamiltonian. 

Such a three-step procedure allows, respectively, for (i) taking into account the 
competition between the nuclear and Coulomb interactions in building up the single 
Slater determinant, which is becoming 'deformed' in the isospace, (ii) restoring the 
isospin symmetry, and (iii) letting the nuclear and Coulomb interactions pick the 
correct mixtures of symmetry-restored eigenstates of the isospin. 

In the present study, we briefly overview the main theoretical building blocks of 
the formalism (Sect.[2|) and discuss preliminary applications. In particular, we give 
results for the isospin-mixing parameters (ac) calculated in the ground states of 
N = Z nuclei (Sect. [3t . Th e ultimate goal will be to perform simultaneous isospin 
and angular momenturrP^ projections and to systematically calculate the isospin- 
symmetry breaking corrections to the Fermi matrix element (fc) for the set of 
nuclei undergoing the superallowed 0+ 0+ Fermi beta decaji^il^. 

2. Theoretical formalism: isospin restoration and Coulomb 
rediagonalization scheme 

The first step and the starting point of our approach is the determination of the 
isospin-symmetry-broken single-particle (s.p.) Slater determinant |HF) calculated 
by using the Hartree-Fock (HF) theory including the isospin-invariant Skyrme {V^) 
and the isospin-symmetry-breaking Coulomb (V^) interactions: 

H^H^ + V'^ where H^^f + V^. (1) 

The isospace-deformed state |HF) admixes higher isospin components T > jT^j: 

|HF)= bT,TAv;T,T,), (2) 

T>\T^\ 

where T and Tz are the total isospin and its third component, respectively, rj labels 
all other quantum numbers pertaining to the |HF) state, and the coefficients bx.r^ 
are such that X]t>|t^| l^r.r^ P = 1- 

In the second step we create the good-isospin states jry; T, Tz) by projecting them 
out from the Slater determinant |HF): 

W,T,Tz)^-^PlTAm- (3) 
In the following, we denote the mixing coefficients 6t.t^ and 

|6t,t. f = imPlr^ |HF) , E^"^^ = (77; T, Tz\H\m T, Tz), (4) 

as being obtained before rediagonalization. In the above formulae, denotes 
the conventional^ SO (3) projection operator reduced to one dimension due to the 
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Tz quantum number conservation, that is: 

PIt^^^^ dp sin (3d^^^^if3)RiP), (5) 

where R{(3) — e^^^'^^ denotes active-rotation operator by the Euler angle /3 in the 
isospace and is the Wigner rf- function^. 

In the third step we mix the projected states, 

\mn,T,)= alrM,T,T,), (6) 

T>|T^| 

and determine the mixing coefficients ^ by diagonahzing Hamiltonian ([T]) in the 
space of projected states, 

J2 (77;T,r,|ff|77;r',T,)aJ,^^^ aj,^^, (7) 

T'>|T^| 

where n enumerates the obtained eigenstates. In the following, we denote the mixing 
coefficients t and eigenenergies E^^_ as being obtained after rediagonalization. 
The lowest-energy solution, for n = 1, corresponds to the isospin mixing in the 
ground state. 

The Skyrme Hamiltonian, , is an isoscalar operator; hence, it contributes 
only to the diagonal matrix elements of the Hamiltonian ([1]), (r/; T, Tz\H^\ri; T,Tz), 
which can be obtained from: 

(HFliJ-^P^^-rjHF) = / d/3sin/34,T,(/3)(HF|iJ^i?(/3)|HF). (8) 
Jo 

Similarly, calculation of the diagonal and non-diagonal matrix elements of the 
Coulomb interaction, {r);T,Tz\V^\T];T' ,Tz), can be efficiently performed after de- 
composing F*^ into the isoscalar, V^, isovector, V^J, and isotensor, compo- 
nents, and by making use of the SO (3) transformation rules for the spherical ten- 

12 

sors under rotations in the isospace . In the particular case of one-dimensional 
projection we deal with in this work, all matrix elements of axial spherical tensors 
reduce to one-dimensional integrals over the Euler angle /3: 

A 

fj,' — — X 

d/3 sin/3 4,^^^(/3) (HF|T/f^,i?(/3)|HF), (9) 

where T!, = — ft' and C'i^Tf, denote standard Clebsch-Gordan coefficients. 
The Skyrme-Hamiltonian and Coulomb-interaction kernels, (HF|_ff'^i?(/3)|HF) and 
(HF|y'-^J?(/3)|HF), respectively, can be evaluated by using expressions for the stan- 
dard diagonal kernelsl^J (/3 = 0) and replacing there the isoscalar and isovector 
densities and currents with the so-called transition densities and currents. Exact 
direct and exchange kern els of the Coulomb interaction can be evaluated by using 
methods outhned in Refs.l^HSHll. 
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Fig. 1. The isospin-mixing parameter ac calculated in ^''"Sn (after rediagonalization ) as a func- 
tion of the number of the HO shells A^o • The results were obtained by using the SIIll^ Skyrme 
parameterization. 




3. Numerical applications: the isospin mixing and the 
isospin-projected energies in N=Z nuclei 

The isospin-mixing parameter, calculated before and after rediagonalization, is de- 
fined as ac = 1 — |^|t^|.tJ^ and ac = 1 — hj^fr 1^' respectively. Its theoretical 
accuracy depends on different factors, and in particular, on the size of the spherical 
harmonic-oscillator (HO) basis used in the calculations. A choice of the number of 
the HO shells Ng included in such calculations is always a result of a trade-off be- 
tween the accuracy and the CPU-time efficiency. In this respect, a bottle-neck in our 
calculation scheme is the exact treatment of the exchange Coulomb contribution, 
which makes calculations prohibitively time consuming for iVo > 16. 

Dependence of the isospin-mixing parameter on A^o is depicted in Fig.[TJ The fig- 
ure shows ac in ^*"^Sn, calculated after rediagonalization, by using the SIII Skyrme 
parameterization of In the expanded scale of the figure, a significant vari- 

ation of the mixing parameter with Nq is clearly seen. Unfortunately, the mixing 
parameter ac does not stabilize at Nq ~ 16. Hence, by using the present method, 
ac cannot be calculated with the absolute precision greater than ±0.002, or with 
the relative precision grater than ±4%. However, our studies show that the inaccu- 
racy in evaluating ac due to the basis cut-off appears to be much smaller than the 
uncertainty related to the Skyrme force parameterizatiorffSI. 

Fig. [2] shows the total binding energies versus Nq, calculated for doubly magic 
nuclei: ^^O, ^"Ca, ^^Ni, and ^"^Sn. This set of calculations was performed by us- 
ing the SII Skyrme force parameterizatiorff^. The curves labeled by black squares 
depict the projected energies E^^^ calculated before rediagonalization, and 
those marked by triangles show the total binding energies 

En=i 0' obtained 

after rediagonalization of the total Hamiltonian in the isospin-projected basis. The 
figure shows that (z) the Coulomb rediagonalization effect increases with increasing 
Z as anticipated, and that {ii) the choice of iVo = 12 HO shells provides a reason- 
able estimate for the total binding energy even for ^'^'^Sn. Hence, all calculations 
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Fig. 2. Total binding energies in ^^O, ^*^Ca, ^^Ni, and ^'^'^Sn as a function s of the number of the 
HO shells A^o- Calculations were performed by using the SII Skyrme iorc^. Dots and squares 
label the binding energies calculated before and after rediagonalization, respectively. 
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Fig. 3. The isospin mixing (left) and the difference between the total binding energy and the HF 
energy (right) in N = Z nuclei, calculated for the SIII Skyrme force and A^o = 12 HO shells. 
The results shown by open and full dots represent variants of the calculation before and after 
rediagonalization, respectively. 



presented below are done for Nq ~ 12. 

Fig. [3K shows the isospin mixing in N = Z nuclei as a function of the mass 
number A. Results obtained before and after rediagonalization are shown by open 
and full dots, respectively. In both variants of the calculations, the isospin mixing 
shows a gradual increase as a function of A. It increases from a fraction of a percent 
in ^^O to about 4%-5%, depending on the variant of the calculation. Note that the 
results obtained before rediagonalization follow closely those obtained in Rem. 

The isospin mixing obtained after removing the spurious mean-field component 
through the Coulomb rediagonalization is systematically larger than the one ob- 
tained within the HF method followed by the exact isospin projection. This result 
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confirms that the mean-field breaks the isospin symmetry in such a way that it coun- 
terbalances the external symmetry breaking mechanism caused by the Coulomb 
field. Nevertheless, as clearly seen in Fig. [Sjj, the HF energy is astonishingly close 
to the total energy obtained after the Coulomb rediagonalization. 

4. Summary 

In the present study, we presented a new theoretical tool that allowed for isospin 
projection of Slater determinants and we discussed first applications of the formal- 
ism to calculate the isospin-mixing parameters ac and total binding energies in 
N = Z nuclei. In particular, we discussed the basis-size dependence of ac and we 
showed that the basis truncation may introduce about ±4% uncertainty in ac- 

We also discussed the role and magnitude of the spurious isospin-symmetry- 
violating response of the self-consistent mean field against the physical symmetry- 
breaking effects of the Coulomb field. We showed that even in N = Z nuclei, 
the self-consistent mean-field may induce unphysical isospin mixing that reduces 
ac by as much as 30% in ^^^Sn. This unphysical mechanism is due to the very 
variational nature of the self-consistent mean-field scheme, which introduces its 
own isospin-symmetry breaking field and partly counterbalances the repulsive 
symmetry-breaking Coulomb field so as to minimize the total binding energy. Nev- 
ertheless, our calculations show that the HF binding energies follow extremely 
closely those obtained by rediagonalizing the Hamiltonian within the set of isospin- 
projected states. 

This work was supported in part by the Polish Ministry of Science under 
Contract No. N N202 328234 and by the Academy of Finland and University of 
Jyvaskyla within the FIDIPRO programme. 
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